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10 Abstract 

1 1 High resolution aerosol, cloud, water vapor, and atmospheric profile data from the Moderate 

12 Resolution Imaging Spectroradiometer (MODIS) are utilized to examine the impact of 

13 aerosols on clouds during the Amazonian biomass burning season in Rondonia, Brazil. It is 

14 found that increasing background column water vapor (CWV) throughout this transition 

15 season between the Amazon dry and wet seasons exerts a strong effect on cloud properties. 

16 Asa result, aerosol-cloud correlations should be stratified by column water vapor to achieve a 

17 more accurate assessment of the effect of aerosols on clouds. Previous studies ignored the 

1 8 systematic changes to meteorological factors during the transition season, leading to possible 

1 9 misinterpretation of their results. Cloud fraction is shown generally to increase with aerosol 

20 optical depth (AOD) for both low and high values of column water vapor, whereas the 

21 relationship between cloud optical depth (COD) and AOD exhibits a different relationship. 

22 COD increases with AOD until AOD ~ 0.25 due to the first indirect (microphysical) effect. 

23 At higher values of AOD, COD is found to decrease with increasing AOD, which may be due 

24 to: (1) the inhibition of cloud development by absorbing aerosols (radiative effect) and/or (2) 

25 a retrieval artifact in which the measured reflectance in the visible is less than expected from a 

26 cloud top either from the darkening of clouds through the addition of carbonaceous biomass 

27 burning aerosols or subpixel dark surface contamination in the measured cloud reflectance. If 

28 (1) is a contributing mechanism, as we suspect, then a linear relationship between the indirect 

29 effect and increasing AOD, assumed in a majority of GCMs, is inaccurate since these models 
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1 do not include treatment of aerosol absorption in and around clouds. The effect of aerosols on 

2 both column water vapor and clouds over varying land surface types is also analyzed. The 

3 study finds that the difference in column water vapor between forest and pasture is not 

4 correlated with aerosol loading, supporting the assumption that temporal variation of column 

5 water vapor is primarily a function of the larger-scale meteorology. However, a difference in 

6 the response of cloud fraction to increasing AOD is observed between forest and pasture. 

7 This suggests that dissimilarities between other meteorological factors, such as atmospheric 

8 stability, are likely to have an impact on aerosol-cloud correlations between different land- 

9 cover types. 

10 

1 1 1 Introduction 

12 The effect of aerosol particles on the hydrological cycle remains one of the largest 

13 uncertainties in our climate system. Biomass burning, from both deforestation and annual 

14 agricultural burning, is the largest anthropogenic source of such particles in the Southern 

15 Hemisphere. A variety of observational and modeling studies have examined the effect of 

16 aerosols on the regional hydrometeorology over the Amazon Basin during the biomass 

17 burning season (Kaufman and Nakajima 1993, Kaufman and Fraser 1997, Koren et al. 2004, 

1 8 Feingold et al. 2005, Yu et al. 2007, Koren et al. 2008, Zhang et al. 2008, Martins et al. 2009). 

19 In these, and in other studies, biomass burning aerosols have been shown to affect clouds 

20 through both microphysical and radiative mechanisms (Kaufman et al. 2005b, Kaufman and 

21 Koren 2006, Koren et al. 2008, Rosenfeld et al. 2008). Depending on the concentration of 

22 aerosol, its chemical composition, size distribution, vertical distribution relative to clouds, and 

23 background cloud characteristics, biomass burning aerosols are suggested either to inhibit or 

24 invigorate cloud formation and/or growth (Feingold et al. 2005, Koren et al. 2008, Chand et 

25 al. 2009). 

26 Carbonaceous biomass burning aerosols can absorb solar radiation, warming the aerosol layer 

27 and reducing the radiation reaching the surface (Ackerman et al. 2000, Kaufman et al. 2002, 

28 Koren et al. 2004). This effect cools the surface, stabilizes the lower troposphere, suppresses 

29 surface heat and moisture fluxes, and slows the hydrological cycle (Jacobson 2002, Andreae 

30 et al. 2004, Ramanathan et al. 2005). Evaporation of clouds within the aerosol layer may also 

3 1 occur due to the increase in temperature and decrease in relative humidity caused by aerosol 

32 absorption of solar radiation (Yu et al. 2002, Jacobson 2006). Referred to as absorptive 
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1 effects or the semi-direct effect, these radiative processes primarily suppress the formation 

2 and growth of clouds. 

3 Microphysical effects, on the other hand, enhance cloud formation and growth. Biomass 

4 burning aerosols are hygroscopic and can serve as cloud condensation nuclei (Feingold et al. 

5 2001, Andreae et al. 2002, Andreae et al. 2004). Expansion chamber experiments have 

6 shown that the addition of cloud condensation nuclei nucleates a larger number of smaller 

7 cloud droplets, and these droplets are therefore slower to coalesce to form precipitation (Gunn 

8 and Phillips 1957, Squires et al. 1958). These aerosol processed clouds are also more 

9 reflective, change drizzle properties, and have been suggested to have longer lifetimes 

10 (Twomey 1977, Albrecht 1989). However, the apparent darkening of clouds by absorptive 

1 1 biomass burning aerosols has also been observed by satellite (Kaufman and Nakajima 1993, 

12 Wilcox et al. 2009). Recent studies have shown that these polluted clouds may become 

13 invigorated, with higher liquid or ice water paths and cloud top pressures (Andreae et al. 

14 2004, Khain et al. 2005, Koren et al. 2005, Lin et al. 2006, Koren et al. 2008, Rosenfeld et al. 

1 5 2008, Meskhidze et al. 2009). The delay of raindrop formation in polluted clouds suppresses 

1 6 downdrafts, which allows for the generation of greater updrafts and stronger convection. The 

17 updrafts carry water vapor to higher altitudes, where additional energy from the latent heat of 

1 8 freezing may be released, further invigorating convection (Williams et al. 2002, Andreae et al. 

19 2004, Rosenfeld et al. 2008). Increases in aerosol optical depth have also been linked to 

20 increases in cloud fraction, particularly at low AODs (Koren et al. 2005, Lin et al. 2006, 

2 1 Myhre et al. 2007, Yu et al. 2007). 

22 More recent studies have illustrated that there may be a smooth transition between these 

23 competing microphysical and radiative effects (Koren et al. 2008, Rosenfeld et al. 2008). 

24 Using MODIS data over the Amazon, Koren et al. (2008) hypothesizes that microphysical 

25 processes dominate at lower AODs, increasing cloud fraction and height, whereas radiative 

26 processes dominate at higher AODs, decreasing cloud fraction and height. The study also 

27 shows that the relative contributions of the microphysical and radiative effects are strongly 

28 tied to the initial cloud fraction - the radiative absorption effect begins to dominate at lower 

29 values of AOD for lower initial cloud fractions. This is due to the hypothesized aerosol 

30 absorption cloud fraction feedback (AFF): Stabilization of the near-surface atmosphere due 

31 to aerosol absorption of radiation initially reduces cloudiness, which then exposes more of the 

32 aerosol layer, further reducing cloudiness (Koren et al. 2008). Eventually, the reduction in 
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1 cloudiness will allow sufficient surface heating to destabilize the lower atmosphere and 

2 reverse the positive feedback (Koren et al. 2008). This reduction in cloudiness has also been 

3 described as a BC-low cloud positive feedback loop (Jacobson 2002). For low cloud 

4 fractions, more of the aerosol layer is available for absorption, resulting in a stronger 

5 feedback. For higher cloud fractions, microphysical invigoration will dominate for the same 

6 degree of aerosol loading. An implication of this finding is that aerosol effects on either 

7 sparse or dense cloud fields may have entirely opposite effects on climate forcing. 

8 For the majority of observational biomass burning studies over the Amazon, the months of 

9 August, September, and October are customarily selected due to the combination of high 

10 aerosol loading from biomass burning and consistent dry conditions present during these 

11 months. The Amazon dry season occurs during the southern hemisphere winter, and is 

12 defined by a subtropical high pressure that moves from the Atlantic Ocean over the Amazon 

13 Basin (Nobre et al. 1998). Felled vegetation is allowed to dry out during the season, and is 

14 burned during the southern hemisphere spring. High pressure typically remains over the 

15 region until mid-to-late October, when the onset of the rainy season begins (Nobre et al. 

16 1998). During this biomass burning transition season, cloud properties are assumed to be 

1 7 weakly dependent on meteorology due to the stationary high pressure overhead (Kaufman and 

18 Nakajima 1993, Koren et al. 2004, Yu et al. 2007, Koren et al. 2008). However, this study 

19 illustrates that gradually increasing background column water vapor (CWV) during the 

20 biomass burning transition season, between the end of the dry season and the beginning of the 

21 wet season, has a discemable impact on aerosol-cloud correlations. 

22 Previous studies have noted that aerosol loading and CWV are weakly spatially correlated 

23 over the Amazon Basin on seasonal timescales (Feingold et al. 2001). On smaller spatial and 

24 temporal scales, latitudinal variation of both CWV and AOD will produce high correlations in 

25 some areas (Kaufman and Fraser 1997). Regions where simultaneous advection of these two 

26 parameters takes place also exhibit high correlations (Remer et al. 1998). However, these 

27 studies do not take into consideration inconsistent sampling resulting from temporal 

28 correlations on longer intraseasonal timescales. Variability in the CWV, which can be 

29 observed from space on high temporal and spatial scales using the MODIS sensor, may be 

30 used as a tracer for large-scale meteorological variability. Attempts to remove the influence 

31 of meteorology from aerosol-cloud correlations are common; a study by Koren et al. (2005) 

32 accounts for meteorology by stratifying data by vertical wind velocity, an NCEP parameter 
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1 that correlates well with MODIS cloud parameters. Here, we use column water vapor, which 

2 is shown also to correlate well with MODIS cloud parameters, to stratify the aerosol and 

3 cloud data. 

4 The latter portion of the paper tests the assumptions made in the former portion of the paper 

5 by analyzing effects of aerosols on clouds and column water vapor over different land surface 

6 types. Many studies have probed the local effects of deforestation on local and regional 

7 meteorology; however no observational studies have analyzed the aerosol effect on clouds 

8 over both forested and deforested land. It is estimated that 50% to 60% of total rainfall in the 

9 Amazon is derived from evapotranspiration (Marques et al. 1977, Salati and Nobre 1991). 

10 Regional modeling studies have estimated a reduction in evapotranspiration of about 30% 

1 1 over deforested regions of the Amazon compared with forested regions (Nobre et al. 1991). 

12 A group of studies has probed into the regional climate effects of deforestation through 

13 changes in surface energy and water vapor fluxes and land-atmosphere interactions 

14 (Henderson-Sellers and Gomitz 1984, Nobre et al. 1991, Cutrim et al. 1995, Wang et al. 2000, 

1 5 Roy et al. 2002, Negri et al. 2004, Wang et al. 2009). Depending on the structure and scale of 

16 the deforestation, contrasting effects on clouds and precipitation are observed (D’ Almeida et 

17 al. 2007). Smaller scale deforestation (i.e. the fish bone pattern) spawns mesoscale 

18 circulations that arise from land surface heterogeneities on the finer scale (Segal et al., 1988, 

19 Wang et al. 2000, Roy et al. 2002, Wang et al. 2009). Enhanced shallow convection over 

20 deforested regions is caused in part by a land breeze from nearby moisture-rich forests. When 

21 this moist land breeze reaches unstable air over the deforested region (due to greater surface 

22 heating), it rises to form clouds (Segal et al. 1988, Roy et al. 2002). Several observational 

23 studies have shown an increase in shallow convection over disturbed regions of the Amazon 

24 due to this direct thermal circulation, particularly in the state of Rondonia, Brazil (Cutrim et 

25 al. 1 995, Negri et al. 2004, Wang et al. 2009). This study differs from previous studies in that 

26 the aerosol effect on water vapor and clouds over forested and deforested is examined, a 

27 logical extension of the studies noted above. 

28 Accounting for meteorological variability in observational aerosol-cloud correlation studies is 

29 paramount. The aerosol effect on cloud properties over the Amazon has been shown to differ 

30 in a humid year compared to a dry year (Yu et al. 2007). In the humid year, cloud fraction 

31 predominantly increases with increasing AOD, whereas in a dry year, cloud fraction 

32 predominantly decreases with AOD. The mechanisms behind this difference are still largely 
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1 unexplained (Yu et al. 2007). This study incorporates both dry and wet years, reducing this 

2 interannual variability, while also stratifying the data by CWV to ensure similar background 

3 moisture conditions exist along the range of AOD retrievals used. This study probes into two 

4 primary questions: (1) how can contamination from meteorological variability in observed 

5 aerosol-cloud relationships be reduced and (2) what is the aerosol effect on clouds and 

6 column water vapor over different land surface types. 

7 

8 2 Data and Methods 

9 MODIS, onboard the Terra and Aqua satellites, provides relatively high spatial resolution 

10 (250 m - 500 m) while achieving near global coverage on a daily basis (Salomonson et al. 

1 1 1989). The Terra satellite is on a sun synchronous orbit with an overpass at approximately 

12 10:30 AM local time, whereas the Aqua satellite overpasses at approximately 1:30 PM local 

1 3 time. The Aqua pass is chosen over the Terra pass since warm clouds are more likely to be 

14 developed in the afternoon compared to the morning. We choose MODIS over other sensors, 

15 such as the Multi-angle Imaging Spectroradiometer (MISR), due to its high resolution and 

1 6 because MODIS produces a variety of cloud and atmospheric profile products, which other 

17 sensors do not (Diner et al. 1998). This study employs MODIS Swath Level 2 aerosol, cloud, 

1 8 water vapor, and stability products (King et al. 2003). Aerosol optical depth from the aerosol 

19 product is calculated over land and ocean at a wavelength of 550 nm, with a footprint of 10 

20 km x 10 km (Kaufman et al. 1997, Remer et al. 2005, Levy et al. 2007). Surface reflectances 

21 utilized in the aerosol algorithm are dynamic, and related to empirical functions which match 

22 reflectances in the near IR to visible wavelengths in order to calculate an NDVI-like measure 

23 of vegetation and geometry (Levy et al. 2007). Validation with ground-based AERONET 

24 observations yield an overall error of ±0.05 ± 0.1 5x a over land, where x a is the aerosol optical 

25 depth at 550 nm (Levy et al. 2007). 

26 The Level 2 cloud product produces retrievals of cloud fraction, cloud top properties, cloud 

27 phase properties, and cloud microphysical properties, calculated using fourteen of the thirty- 

28 six MODIS spectral bands (Ackerman et al. 1998). The cloud fraction and cloud top 

29 properties are produced at 5 km x 5 km resolution, whereas the microphysical properties are 

30 produced at 1 km x 1 km resolution (Platnick et al. 2003). The 5 km cloud fraction product is 

3 1 computed using the fraction of 1 km cloudy pixels in the 5 km footprint, as determined by the 

32 cloud mask product (Platnick et al. 2003). The cloud mask product is based on the principle 
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1 that clouds are characterized by higher reflectance and lower temperature than the earth’s 

2 surface, and is clear-sky conservative (Ackerman et al. 1997). The cloud phase retrieval is 

3 determined by analyzing the contrasting effects of water droplets and ice crystals on the 

4 brightness temperature in the infrared bands (Platnick et al. 2003). 

5 The cloud microphysical properties include cloud optical depth (thickness), cloud effective 

6 radius, and cloud liquid water path. Cloud optical depth (COD) is inversely calculated from 

7 spectral reflectance measurements and surface albedo data by using a lookup table (Nakajima 

8 and King 1990, Twomey and Cocks 1982, Twomey and Cocks 1989). For over-land 

9 retrievals, COD is calculated at wavelengths 0.645 pm and 2.130 pm. The 1 km and 5 km 

10 cloud data are averaged into 10 km x 10 km grid boxes in order to conform to the Level 2 

1 1 aerosol data. 

12 Column water vapor (precipitable water vapor) is derived by integrating the 101 levels at 

13 which water vapor mixing ratios are calculated in the MOD07 atmospheric profile product 

14 (Seemann et al. 2002). Units are reported in centimeters of equivalent liquid water. We 

1 5 choose not to use the near-IR column water vapor product due to its limitations over dark 

16 surfaces (Gao and Kaufman 1998). The profile algorithms are based on thermal emissions of 

17 atmospheric gases with uniform distributions, such as oxygen and carbon dioxide (Seemann 

18 et al. 2002). The 1000 hPa temperature is calculated using the skin temperature, surface 

19 pressure, and Poisson’s Equation for potential temperature. Skin temperature is not the best 

20 approximation for surface air temperature; however, for the purposes of finding relative 

21 differences in atmospheric stability across our study domain, we feel this is our best option 

22 considering that a high-resolution observational temperature data set is currently not available 

23 from any other source. The water vapor and profile products are also produced at 5 km 

24 resolution, but are averaged to the 10 km scale to conform to the Level 2 aerosol data. 

25 Products requiring a clear sky, such as the MOD07 products, are able to be estimated for a 

26 majority of 1 0 km pixels with cloud fractions less than one. 

27 In the first portion of the paper, aerosol and cloud data over all non-water surfaces are 

28 included in the analyses. In the latter portion of the paper, MODIS aerosol, cloud, stability, 

29 and water vapor data are stratified by land cover type. This requires an up-to-date, high- 

30 resolution land cover classification data set. The Land Processes Distributed Active Archive 

31 Center (LP DAAC), located at the U.S. Geological Survey (USGS) Earth Resources 

32 Observation and Science (EROS) Center, provides a combined Terra/Aqua yearly land cover 
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1 product - MCD12Q1 (http://lpdaac.usgs.gov). This product employs MODIS BRDF-adjusted 

2 surface reflectances, land surface temperature data, enhanced vegetation index data, and 

3 terrain elevation information along with neural network classification algorithms and training 

4 data to assign land cover classifications (Strahler et al. 1999). The data are resampled from 

5 500 m x 500 m resolution to 0.1° x 0.1° resolution to approximately match the resolution of 

6 the Level 2 swath aerosol, cloud and profile data. 

7 Our 5° x 5° study region encompasses the deforested region of Ji Parana in Rondonia, Brazil, 

8 as well as a protected forest to the east, illustrated in Fig 1 . Green regions represent primary 

9 forest and yellow regions represent pasture. Blue regions represent land cover classes that do 

10 not explicitly fall in either class. Broad-leaf forest classifications are assigned to the forested 

1 1 category, whereas closed shrublands, open shrublands, woody savannas, savannas, grasslands, 

12 croplands, cropland and natural vegetation mosaic, and barren or sparsely vegetated 

1 3 classifications are assigned to the pasture category, according to the International Geosphere- 

14 Biosphere Programme categorization scheme (Strahler et al. 1999). The percentage of pasture 

1 5 increases with time in our fixed study region due to ongoing deforestation. 

16 A small region is chosen for this study, compared with other studies of its type, so that 

17 meteorological differences due to spatial variation will be better removed, and so that 

18 stratification of atmospheric data by land cover type can be conducted. The high spatial 

19 resolution of the Level 2 data allows for the accumulation of a sufficient data record for 

20 analysis, while also reducing the pixel size in which aerosol and cloud characteristics need to 

21 be assumed constant. Previous studies have relied on Level 3 data, which assumes aerosol 

22 and cloud characteristics are constant over 1° x 1° pixels (Koren et al. 2008). The study 

23 encompasses the Amazonian biomass burning months of August, September and October for 

24 the years 2004 through 2007. This time period includes both dry and wet years, as well as 

25 years with both heavy and light biomass burning, according to MODIS and NCEP/NCAR 

26 Reanalysis data (Kalnay et al. 1996, Koren et al. 2007). The inclusion of multiple years of 

27 data reduces the interannual variability present in the data. In addition, using multiple years 

28 allows for the compilation of a large enough data set for analysis in the small study region. 

29 To account further for meteorological variability, NCEP/NCAR Reanalysis 700 hPa wind 

30 vectors are used to remove days for which the South Atlantic Subtropical High was not the 

31 dominant weather pattern over the region (Kalnay et al. 1996). This resulted in a removal of 

32 43 days over the four year period, equaling 12% of the total number of days in the months of 
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1 August through October in 2004 through 2007. All retrievals that are not considered “useful” 

2 or were considered “bad” quality according to the Level 2 quality assurance bit data were also 

3 removed. Warm clouds are segregated from cold and unknown-phase clouds by only 

4 retaining 10 km pixels that contain >95% of 1 km cloud retrievals in the liquid water phase. 

5 The re-sampled 10 km atmospheric aerosol, cloud, and profile data are then segregated by 

6 forest and pasture in the latter portion of the paper to explore the effect of land-atmosphere- 

7 aerosol interactions on CWV and clouds. 

8 Aerosol optical depth and column water vapor MODIS satellite retrievals are compared with 

9 Aerosol Robotic Network (AERONET) sun photometer data in Fig SI (Holben et al. 1998). 

10 Data represent two stations within the 5° x 5° study region, Abracos Hill, active in 2004 and 

1 1 2005, and Ji Parana SE, active in 2006 and 2007. Only days used in the study are included. 

12 Correlation coefficients between MODIS and AERONET retrievals are above 0.90 for both 

1 3 aerosol optical depth and column water vapor, providing confidence in the MODIS satellite 

14 retrievals. 

15 

16 3 Results 

17 3.1 Effect of Intraseasonal Meteorological Variability 

1 8 CWV increases over our study region throughout the dry-to-wet transition season between the 

1 9 months of August to October. Fig 2a shows MODIS CWV and liquid water cloud fraction 

20 averaged into eight equally-spaced bins between August 1 st and October 31 st , for all years 

21 between and including 2004 and 2007. Cloud fraction and CWV are understandably highly 

22 positively correlated, as water vapor is one of the principal components required for cloud 

23 formation. Fig 2b shows MODIS aerosol optical depth at 550 nm, also binned by day of the 

24 season, similar to Fig 2a. Unlike CWV, AOD increases from August 1 st until approximately 

25 the middle of September, at which point the trend reverses and decreases until the end of the 

26 season on October 31 st . This mid-season peak in aerosol loading is indicative of a biomass 

27 burning peak that is largely determined by social behavior. Burning customarily occurs late 

28 enough into the season to sufficiently allow vegetation to dry out, but not too late to risk an 

29 early onset of the rainy season. 

30 Fig 2 indicates that CWV and cloud fraction are directly correlated with aerosol optical depth 

3 1 during the first half of the season, whereas these variables are inversely correlated during the 
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1 second half of the season. Thus, if we plot CWV and cloud fraction directly against AOD, 

2 separated into the two halves of the season, we would expect the signs of the regression 

3 slopes to be the same for each parameter in each half season. If CWV increases with aerosol 

4 in the first half season, then so should cloud fraction. If CWV and AOD are negatively 

5 correlated in the second half season, then cloud fraction should also be negatively correlated 

6 with AOD. 

7 Fig 3a shows CWV binned by aerosol optical depth for all non-zero cloud fraction retrievals 

8 in our study region for the same time period as in Fig 2. The data are stratified by day of the 

9 season, with plots representing each half of the season and all days in the season. Fig 3b 

10 depicts liquid cloud fraction binned by aerosol optical depth, for the same periods as in Fig 

11 3a. In these figures, cloud properties and CWV are binned by AOD, with each bin 

12 representing 12.5 percentile of the AOD values. This method has been used previously in 

1 3 other studies so that bias is not introduced through inconsistent sampling in each bin (Lin et 

14 al. 2006). The absolute low and high AOD boundaries are assigned to be 0.05 and 0.8, 

1 5 however, the location of individual bin edges vary for each plot. The bin centers are defined 

16 as the average AOD value in each bin. Error bars representing the standard errors of the bin 

17 average (o/N ) are included for each bin. Because the number of samples is equal for each 

18 bin due to the sampling procedure used, the standard error is directly proportional to the 

1 9 standard deviation of the samples in each bin. 

20 Even though AOD values greater than 0.8 are routinely observed in this region during the 

21 biomass burning season, higher values are not included to prevent aerosol misclassification as 

22 cloud (Brennan et al. 2005). In this range, significant cloud contamination of aerosol by 

23 clouds is also unlikely, yet occasional contamination may still occur (Kaufman et al. 2005b). 

24 It is also improbable that these correlations are due to a 3D cloud effect, which artificially 

25 increases AOD retrievals in the regions neighboring clouds (Wen et al. 2006). This effect has 

26 been shown to be larger for greater cloud cover and aerosol loading, but our results indicate 

27 that the strongest positive correlation between cloud fraction and AOD occurs at low AODs 

28 (Yu et al. 2007). It is also unlikely that the positive cloud ffaction-AOD correlation is due to 

29 atmospheric stability, since increasing AOD has been shown to correlate with increasing 

30 stability (Davidi et al. 2009). 

31 In the first half season, Aug. 1 st to Sep. 1 5 th , there is a positive correlation between both CWV 

32 and AOD, as well as between cloud fraction and AOD, as expected from Fig 2. However, in 


10 



1 the second half season, Sept. 15 th to Oct. 31, the CWV is negatively correlated with AOD as 

2 expected, but cloud fraction is not. This deviation of the cloud fraction relationship to AOD 

3 in a direct correlation overrides the controlling effect of the seasonal CWV evolution on cloud 

4 formation and strongly points to aerosol as a modifying element of cloud properties. We 

5 conclude that both large-scale meteorological factors traced by seasonally varying CWV and 

6 aerosol loading contribute to changes in water cloud fraction in our region of interest. 

7 This conclusion has important implications for aerosol-cloud analysis. Following previous 

8 studies, common practice would be to accumulate all pixels with both aerosol and cloud 

9 retrievals during the August-October season, and then bin the cloud retrievals by aerosol 

10 optical depth, implicitly assuming constant meteorological conditions throughout this period 

11 (Koren et al. 2004, Yu et al. 2007, Koren et al. 2008). However, this study finds that 

12 embedded in these correlations are systematic variations in CWV with aerosol optical depth. 

13 As a result, this intraseasonal CWV signal may impact aerosol-cloud correlations. For 

14 example, the “boomerang” shape identified by Koren et al. (2008) and attributed to the 

1 5 combination of microphysical and radiative effects by aerosols on clouds may also contain an 

1 6 element of intraseasonal evolution of the meteorological conditions. In Fig 3a we identify a 

17 boomerang shape in the CWV vs. AOD plot that is due solely to the combination of the 

1 8 different halves of the season, and not a physical consequence of the aerosol at all. The 

1 9 similar boomerang shape in cloud fraction over the full biomass season seen in Fig 3b reflects 

20 both the actual relationship between cloud fraction and aerosol in the latter half season and 

21 also the artifact created by combining the two half seasons. 

22 Fig 3c and Fig 3d show CWV and cloud fraction again binned by AOD, but only for cloud 

23 fractions less than 0.5. By retaining only low cloud fractions, more of the aerosol layer is 

24 exposed to sunlight in each 10 km x 10 km retrieval box. In these scenes, it has been 

25 hypothesized that the increased aerosol exposure to sunlight will result in a stronger 

26 absorption effect on clouds and that the decrease of cloud fraction with AOD will begin at 

27 lower values of aerosol loading (Koren et al. 2008). The lower cloud fraction should not 

28 affect the relationship between CWV and AOD, because CWV should not be affected by a 

29 warming atmospheric column. Indeed, Fig 3c shows little difference in both the shape of the 

30 graph and the magnitude of the change of CWV with AOD when compared with Fig 3 a. 

3 1 Exposure to sunlight has little effect on column water vapor. However, consistent with Koren 

32 et al. 2008, we find that exposure to sunlight does affect the relationship between aerosol 


11 



1 loading and cloud fraction. The less cloudy plot of Fig 3d has shifted the turning point 

2 between increasing and decreasing cloud fraction to lower aerosol loading. Fig 3d locates the 

3 turning point at AOD = 0.30, as opposed to the total data set of Fig 3b where the turning point 

4 is at AOD = 0.55. 

5 Fig 2 and Fig 3 both confirm the associations between clouds and aerosols identified from 

6 satellite retrievals seen in previous studies and also sound a warning that some of these 

7 associations contain an artifact introduced from slowly evolving meteorological factors that 

8 can be traced by CWV. We note a difference between this study and the previous studies 

9 mentioned above. We are focusing on a very small (5° by 5°) box in Ronddnia, and previous 

10 studies took a broader view that included the entire Amazon Basin with its varied surface 

1 1 types, biomes and climatic zones. The trends shown in Fig 2 are applicable for larger areas as 

12 well, but the range between low and high CWV and low and high AOD over the season 

1 3 becomes diluted as the study area is expanded. Our study area captures the strong seasonal 

14 variation in water vapor, cloud and aerosol properties that still exists but in diluted 

1 5 magnitudes when our region is combined with the surrounding areas. (So, as the study area 

1 6 is expanded, the range of CWV and AOD values decreases? This is confusing to me) 

17 We can identify intraseasonal meteorological variations by using CWV as a tracer for the 

18 slow onset of the rainy season and stratify the data set using this parameter to control for 

1 9 meteorology. Fig 4a depicts CWV binned by AOD for all cloud fraction retrievals, and for 

20 five percentile groupings of CWV. Each grouping spans 16 percentile points, with a 

21 minimum percentile of 10% and a maximum percentile of 90% to avoid outlier values. These 

22 bounding low and high percentiles refer to CWV values of 2. 1 0 cm and 5.33 cm, respectively. 

23 In general, the lower percentiles represent retrievals early in the season whereas higher 

24 percentiles represent retrievals late in the season. In each grouping, CWV only varies 

25 marginally between AOD bins, with a maximum difference of 0.09 cm between any two bins 

26 in any grouping. This nominal variation in CWV within each grouping is likely to have only 

27 a minimal effect on cloud properties. 

28 Fig 4b depicts cloud fraction binned by AOD for the five CWV groupings. The change in 

29 cloud fraction with AOD in Fig 4b is more representative of the true aerosol effect on cloud 

30 fraction compared to Fig 3 since the influence of varying CWV has been minimized. For 

31 each CWV grouping, cloud fraction either exhibits a modest increasing trend, on average, 

32 across the AOD range between 0.05 and 0.8, or remains relatively constant along this range. 
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1 The 42 nd to 58 th percentile grouping and the 58 th to 74 th percentile groupings do not exhibit an 

2 increasing trend to the extent of the other groupings. The difference in the aerosol effect on 

3 cloud fraction does not appear to be largely dependent on CWV, which agrees with previous 

4 studies (Feingold et al. 2001). 

5 In general, cloud fraction appears to increase with AOD at the greatest rate at low AODs, 

6 specifically for the lowest CWV grouping as well as the highest CWV grouping. The highest 

7 CWV group, the 74 th to 90 th percentile, demonstrates the greatest microphysical effect - cloud 

8 fraction increases at the greatest rate below an AOD of 0.35 compared to other groupings. In 

9 high water vapor loading environments with low cloud condensation nuclei, addition of 

10 aerosol may increase cloud fraction to a greater extent than in lower water vapor 

11 loading environments due to the impact of stronger updrafts on the lifetime effect 

12 (Albrecht 1989, Kaufman and Fraser 1997) (Can I say this?). This result also agrees with 

13 the aerosol absorption cloud fraction feedback - the greater the initial cloud fraction, the 

14 greater the microphysical effect (Koren et al. 2008). The flattening of these curves at higher 

1 5 AODs suggests a saturation of the microphysical effect. The absence of an aerosol absorption 

1 6 effect in Fig 4b conflicts with some studies (Koren et al. 2004, Koren et al. 2008), but agrees 

1 7 with others (Lin et al. 2006, Yu et al. 2007). Part of the discrepancy between these results and 

1 8 previous ones may be due to the gradient of aerosol absorption properties north to south. The 

19 focus area of this study is embedded in the deforestation zone, with smoke having higher 

20 single scattering albedo (less absorption) than the smoke in the Cerrado to the south (Dubovik 

21 et al. 2002). A less absorbing smoke will have less radiative effect than the model formulated 

22 by Koren et al. (2008). 

23 If the aerosol is affecting cloud microphysics, the signal should be apparent in the cloud 

24 optical depth (COD) as well in the cloud fraction. Fig 5a shows COD binned by AOD for 

25 four out of the five CWV groupings in Fig 4. The lowest CWV grouping is not included due 

26 to relatively large standard errors in each bin resulting from an insufficient sample size. 

27 According to the aerosol first indirect effect, increases in aerosol number loading increase the 

28 number concentration of cloud condensation nuclei, which in turn increase cloud reflectivity 

29 (Twomey 1977). This effect would be observed as an increase in COD with increasing AOD 

30 in Fig 5a. However, Fig 5a illustrates that for all CWV groupings, COD increases only to a 

31 certain AOD threshold between 0.2 and 0.3, and then decreases with increasing AOD. The 

32 highest magnitude peak in COD occurs for the lowest CWV grouping, which also 
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1 corresponds to the lowest cloud fraction from Fig 5. The lowest COD and least negative 

2 correlation with AOD corresponds to the highest CWV and largest cloud fraction. Note that 

3 Fig 5a includes clouds at all stages of vertical development, as long as the clouds are in the 

4 liquid phase. 

5 Fig 5b also shows COD binned by AOD, but only for pixels with average cloud top pressures 

6 between 800 hPa and 850 hPa. This range is roughly between the median and mean liquid 

7 water cloud top pressures over the four year period. Fig 5b shows a similar boomerang 

8 pattern of COD versus AOD as Fig 5a, except with more noise. Only a narrow range of cloud 

9 top pressures are retained to ensure that the COD versus AOD relationships observed in Fig 

10 5 a are not merely the result of sampling clouds at different stages of growth. Even if this 

1 1 were the case, the addition of biomass burning aerosols has been shown to increase cloud top 

12 pressure, which would in turn likely increase COD with increasing AOD, not decrease COD 

13 as shown in Fig 5a (Koren et al. 2005, Koren et al. 2008). The AOD turning point of the 

14 boomerang occurs at generally higher AODs than in Fig 5a, specifically for the lower two 

1 5 CWV groupings. The magnitude of the decrease in COD with increasing AOD in terms of 

16 percentage reduction from the peak is more similar among CWV groupings than in Fig 5a. 

17 This is because in Fig 5b, all clouds have similar cloud top pressures and are situated in more 

1 8 homogenous cloud fields. 

1 9 The increase in COD with AOD for low aerosol loading can be explained by Twomey (1 977) 

20 and follows from aerosol particles affecting the microphysics of the clouds. The decrease in 

21 COD with increasing AOD above AODs of 0.3 in Fig 5a can be explained by a combination 

22 of physical processes and satellite retrieval artifacts. Physically, aerosol absorption of solar 

23 radiation can evaporate or thin the clouds optically. However, satellite retrieval artifacts can 

24 also play a role. Dark smoke above or within clouds (Kaufman and Nakajima, 1993), or 

25 subpixel holes in the clouds that reveal a dark surface beneath reduces the visible reflectance 

26 received by the satellite, and this is interpreted by the retrieval as a lower COD. In fact a 

27 physical-artificial feedback can be started in which the aerosol begins to thin the cloud via 

28 radiative processes, revealing more holes that introduce darker visible reflectance to the 

29 satellite measurements. The microphysical effect dominates at lower AODs whereas the 

30 physical and artificial effects that will decrease COD with aerosol loading are radiative in 

31 nature and will dominate at higher AODs. Considering that all explanations for COD 

32 decrease are dependent on aerosol radiative absorption, the strength of the COD decrease 
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1 should be inversely proportional to cloud fraction, because lower cloud fraction allows the 

2 aerosol greater exposure to sunlight. When the aerosol is exposed to sunlight, more aerosol 

3 translates to even higher rates of absorption and heating. The highest magnitude COD 

4 decrease occurs for the lowest CWV grouping with the lowest cloud fraction, and the lowest 

5 magnitude COD decrease occurs for the highest CWV grouping with the highest cloud 

6 fraction. Furthermore, the AOD threshold where COD changes from increasing to decreasing 

7 occurs at the lowest AOD value for the lowest CWV grouping, and the highest AOD value for 

8 the highest CWV grouping. These two patterns are consistent with the hypothesis that the 

9 COD decrease is a radiative effect of the aerosols, whether the effect is physical or artificial or 

10 a combination of both. 

1 1 Sorting out the exact cause of the COD decrease with AOD is difficult. The absence of a 

12 strong absorption effect in the cloud fraction plots would suggest that darkening may play a 

13 role in the COD plots. Wilcox et al. (2009) found that biomass burning aerosols over a 

14 stratocumulus deck will artificially reduce observed cloud optical depths by the MODIS 

1 5 sensor (Wilcox et al. 2009). The darkening bias was found to be greater for higher AODs and 

16 also for higher non-polluted values of cloud optical depth. Both of these relationships are 

17 observed in Fig 5a. Unlike the Wilcox et al. (2009) study, the data in Fig 5 do not consist 

1 8 only of scenarios with aerosols above clouds. Analysis of individual days using both MODIS 

1 9 Aqua and the CALIPSO lidar, which has an overpass time roughly 1 minute after the Aqua 

20 satellite, results a similar boomerang trend in COD versus AOD even when the cloud layer 

21 forms at the top of the boundary layer, above the aerosols. Fig S2 shows the results of one of 

22 these days, August 12 th , 2006, where clouds form above the aerosol layer and COD exhibits a 

23 boomerang trend with AOD. Colors in the Fig S2b scatter plot represent the relative 

24 uncertainty in each COD retrieval. Low COD uncertainties at high AODs suggest that the 

25 reversal in the COD trend may not be explained solely by uncertainties in the COD retrieval. 

26 Kaufman and Nakajima (1993) also hypothesized that biomass burning aerosol darkens 

27 Amazonian clouds due to the presence of black carbon inside clouds (Kaufman and Nakajima 

28 1 993). This study found that cloud reflectance at 640 nm is reduced from 0.71 to 0.68 for an 

29 increase in aerosol optical depth between 0.1 and 2.0. This small decrease in visible 

30 reflectance (4%) is likely not the sole cause of decreases in COT between 20% and 50%, 

3 1 observed in Fig 5a. Because of the small role expected by the darkening of the clouds and the 

32 absence of a strong absorption effect in the cloud fraction plots, the physical-artificial 
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1 feedback involving thinning clouds and increasing inclusion of dark reflectance from the 

2 surface beneath is the most promising explanation. Modeling studies have also found that 

3 black carbon absorption in clouds results in a non-negligible feedback to surface and mid- 

4 tropospheric temperatures globally (Jacobson 2006). This in-cloud absorption may constitute 

5 up to 10% or more of the total temperature feedback from black carbon. Thus, aerosol 

6 absorption in clouds may have a large effect on the radiative balance of the atmosphere in the 

7 Amazon region. 

8 3.2 Effect of Land-Atmosphere Interactions 

9 Biomass burning in the Amazon affects clouds in two ways - through smoke effects and 

10 through land use changes. Smoke particles can change cloud microphysics and heat the 

1 1 atmospheric cloud environment by absorbing sunlight, as described above. Land use changes 

12 that convert forest to pasture and cropland change the surface heat fluxes, which also affect 

1 3 cloud development. By dividing our data set into forested and deforested (pasture) subsets, 

14 we can begin to understand the relative effects of aerosols and land use changes on cloud 

1 5 properties, and especially how land use changes modify aerosol effects. 

16 Table 1 shows average values of cloud, aerosol, and profile parameters retrieved from the 

17 Level 2 MODIS products. Water cloud fraction is substantially higher for the pasture 

18 compared to the forest, which agrees with previous studies (Cutrim et al. 1995, Negri et al. 

1 9 2004, Wang et al. 2009). In addition, a lower average cloud top pressure is observed for the 

20 pasture, which suggests more shallow cloud development over the pasture compared to the 

21 forest. This result also agrees with previous observational and modeling studies (Wang et al. 

22 2000, Chagnon et al. 2004, Correia et al. 2007, Wang et al. 2009). A consequence of deeper 

23 warm clouds over the forest is larger average liquid water paths over the forest. Column 

24 water vapor is also higher over the forest, due to the reduction in evapotranspiration over the 

25 pasture compared to the forest (Salati and Nobre 1991). Higher atmospheric low-level 

26 stability is observed over the pasture, a parameter defined as the temperature at 850 hPa 

27 minus the temperature at 1000 hPa. The shorter roughness length and lower specific heat of 

28 the pasture result in greater surface heating, as supported by the higher 1000 hPa temperature 

29 compared to the forest. This heating is hypothesized to help spawn shallow cumulus clouds 

30 (Negri et al. 2004). The difference in low-level atmospheric stability, roughly 4 K, correlates 

31 well with other studies (Polcher and Laval 1994, Correia et al. 2007). Aerosol optical depth is 

32 almost identical between the two land cover types, as we would expect since aerosol 
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1 concentrations are often relatively spatially homogenous for hundreds of kilometers, 

2 particularly at some distance away from concentrated smoke plumes (Artaxo et al. 1998, 

3 Smirnov et al. 2000). 

4 Comparison of aerosol effects on CWV and clouds between different land cover types 

5 requires careful consideration of consistent sampling procedures. To perform this analysis, 

6 the same number of samples must be retrieved between forest and pasture for a given day for 

7 each AOD bin to ensure that weighting of the data by time of the year is consistent between 

8 both land cover types. The sampling process is as follows: First, for each parameter and each 

9 AOD bin, the number of samples for each day is determined for forest and pasture. For the 

10 land cover type (forested or deforested) with the minimum number of samples, the same 

1 1 number of samples is taken from the other land cover type. Because bin edges must be 

12 assigned beforehand using this procedure, the number of samples in each AOD bin is not 

13 consistent along all bins. The eight bins for both forest and pasture are equally spaced 

14 between 0.05 and 0.8. 

15 Several randomized mixed forest and pasture simulations are also conducted to assess the 

1 6 significance of the segregated land cover analysis. In each random simulation, the retrievals 

17 are randomly sampled between forest and pasture according to the relative proportion of 

1 8 forest and pasture retrievals in the actual data, so in effect there is no systematic segregation 

19 by land cover in each randomized simulation. Again, the data is processed such that the same 

20 number of samples is retained between forest and pasture for a given day for each AOD bin. 

21 The black line in Fig 6a illustrates the difference in the average sampled Julian Day between 

22 pasture and forest in each AOD bin. The difference is zero for all bins, indicating that the 

23 sampling procedure by definition does not weight the result by day of the year, thus there is 

24 no difference in large-scale meteorology between the pasture and forest for each AOD bin. 

25 The red dotted lines represent the results of five randomized simulations. Again, because the 

26 algorithm does not weight the result by day of the year, there is no difference in large-scale 

27 meteorology between the forest and pasture for each AOD bin in the randomized simulations. 

28 One of the primary underlying assumptions in this study is that aerosols do not impact CWV 

29 through changes in evaporation and transpiration, but instead that CWV is an indicator of the 

30 synoptic scale meteorology. Aerosols could, in theory, affect evapotranspiration rates by 

31 changing both the magnitude and diffuse fractionation of radiation reaching the surface 

32 (Jacobson 2002, Roderick et al. 2001, Min 2005, Knohl et al. 2008, Still et al. 2009). Because 
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1 deforested land has a leaf area index several times lower than forested land, in addition to a 

2 significantly lower heat capacity, aerosols could affect the evapotranspiration and evaporation 

3 from the ground differently between forested and deforested land. We test this assumption by 

4 analyzing aerosol-CWV relationships separately over the forested area and deforested area in 

5 Fig 1. Table 1 indicates that the forest CWV is on average 0.2 cm higher than the pasture 

6 during the dry seasons of 2004 to 2007, which we assume is the result of differences in 

7 evapotranspiration. If we observe this difference correlated with aerosol, then we know that 

8 aerosol is modifying the CWV values and our assumption of using CWV as a tracer for the 

9 large-scale meteorology is incorrect. 

10 Fig 6b shows CWV binned by AOD over the pasture minus CWV binned by AOD over the 

1 1 forest, again for retrievals that also have cloud fractions above zero. Error bars represent the 

12 square root of the sum of the pasture and forest standard errors squared. Red dotted lines 

1 3 again represent the results of the random simulations. There is a lack of a significant trend in 

14 CWV differences with AOD. Instead, CWV is roughly 0.1 cm greater over the forest 

15 compared to the deforested land for all bins, an offset owed to the difference in 

1 6 evapotranspiration between pasture and forest. The random simulations show no such offset 

17 since forest and pasture retrievals are randomly mixed in these simulations. In addition, the 

1 8 magnitude of the variability between AOD bins is similar between the data and the random 

1 9 simulations, further supporting the absence of correlation between CWV and AOD. 

20 Fig 6b shows that addition of aerosols below an AOD of 0.8 does not influence 

21 evapotranspiration sufficiently to affect CWV differently between pasture and forest. 

22 Difference values between pasture and forest do not vary greatly along AOD bins, and the 

23 little variation that does exist is well within the standard error. More generally, the addition 

24 of aerosol below an AOD of 0.8 does not have a noticeable impact on CWV. Thus, changes 

25 in CWV may be attributed completely to changes in synoptic-scale conditions. 

26 Fig 6c shows the difference between cloud fraction binned by AOD over the pasture and over 

27 the forest. Absolute values of aerosol-cloud correlations over forest and pasture are not 

28 computed as they will be influenced by variations in CWV, similar to Fig 3. However, 

29 because the same number of retrievals is retained over the forest and pasture for each day in 

30 each AOD bin, computing cloud fraction differences between forest and pasture in each AOD 

3 1 bin will remove influences of the larger-scale meteorology. 
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1 The cloud fraction difference between pasture and forest is always positive in Fig 6c, 

2 indicating that cloud fraction is greater over the pasture than the forest, as suggested by Table 

3 1 . The figure also illustrates that there are increasingly more clouds over the pasture with 

4 increasing AOD, up to an AOD of 0.5. At this AOD threshold, the trend begins to flatten out, 

5 and even reverses slightly. The five random simulations do not show higher average cloud 

6 fractions over the pasture, nor do they show a stronger microphysical aerosol forcing over the 

7 pasture, supporting the significance of these results in the actual data. What causes this 

8 noticeable dissimilarity in the aerosol effect between pasture and forest at low AODs below 

9 0.5? Because differences in CWV are likely not a driving factor due the results of Fig 6b, it is 

10 suggested that differences in low-level atmospheric stability may be at the root of the 

1 1 dissimilarity. 

12 Fig 6d again shows the difference in cloud fraction binned by AOD between the pasture and 

1 3 forest (black line), but also shows lines representing cloud fraction difference stratified by 

14 stability. The red line includes all retrievals where the low-level stability, defined as the 

15 temperature at 850 hPa minus the temperature at 1000 hPa, is less than -16 K over the pasture 

16 and greater than -16 K over the forest. The blue line includes all retrievals where the low- 

17 level stability is greater than -16 K over the pasture and less than -16 K over the forest. This 

18 threshold stability value of -16 K is selected because it is very close to the median stability 

1 9 over both pasture and forest. In addition, a temperature difference of 1 6 K or greater between 

20 850 and 1000 hPa conservatively corresponds to an unstable atmosphere, according to a 

21 tropical standard atmosphere (A height of 1529 meters corresponds to 850 hPa at 15° N). 

22 The red line represents the difference in cloud fraction between pasture and forest only for 

23 retrievals where the pasture is considered absolutely unstable and the forest is relatively more 

24 stable. In this case (red line), the difference in cloud fraction between pasture and forest 

25 responds similarly to increasing AOD as the case that includes all stabilities (black line). This 

26 similarity is expected since the pasture is often more unstable than the forest. The blue line 

27 represents the reverse situation, retrievals where the forest is considered absolutely unstable 

28 and the pasture relatively more stable. Because this scenario is less likely, larger errors are 

29 observed. Unlike the black and red lines that show increasing clouds over the pasture with 

30 increasing AOD, the blue line does not exhibit such a trend. The absence of a trend may be 

3 1 attributed to the higher stability over the pasture compared to the forest - the microphysical 

32 forcing of clouds may be weaker in atmospheres that lack sufficient low-level instability. 
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1 Low level stability appears to be a factor or a tracer for how aerosol modifies cloud fraction 

2 differently over forested and deforested surfaces. The greater instability that generally occurs 

3 over pasture appears to encourage increases in cloud fraction with AOD. Yet, the response of 

4 COD to increasing AOD, when analyzed, did not appear to be substantially different between 

5 forest and pasture. The statistical approach used here is extremely limited, and cannot 

6 characterize the full mesoscale circulations thought to be responsible for the general 

7 differences in cloud fraction between forested and deforested surfaces (Wang et al. 2009). A 

8 proper 3D cloud resolving model with adequate similations of surface fluxes and mesoscale 

9 circulations will be needed to explain the associations identified in Fig 6c and Fig 6d. 

10 

1 1 4 Conclusions 

12 In this study, a 5° NE x 5° WE region over Rondonia, Brazil was analyzed using high 

13 resolution aerosol, cloud, water vapor, and atmospheric profile data from the Moderate 

14 Resolution Imaging Spectroradiometer (MODIS). Four years of data (2004 - 2007) during the 

15 biomass burning transition season (Aug - Oct) were compiled to analyze the effect of aerosols 

16 on warm cloud development. Several years of data were employed to gather a large enough 

17 dataset for analysis, and to smooth out interannual variability. The MODIS water vapor and 

1 8 aerosol products illustrate that column water vapor (CWV) increased throughout the biomass 

19 burning season, while aerosol loading peaked during the middle of the season. These 

20 intraseasonal trends may produce false correlations in cloud parameter versus AOD plots 

21 when all data throughout the biomass burning season are accumulated and analyzed together. 

22 By dividing the total period into subperiods, we found evidence that both large-scale 

23 meteorological factors traced by CWV and aerosol loading contributed to modification of the 

24 cloud fraction in our area of interest. To better account for this meteorological variability, 

25 data were stratified by column water vapor. When background CWV variability was 

26 removed, cloud fraction increased with AOD, at nearly all levels of CWV, but the strongest 

27 increases occurred at the highest values of column water vapor. 

28 Decrease of cloud fraction with AOD, attributed to aerosol absorption effects, was not 

29 observed in the cloud fraction versus AOD plots once the data were stratified by CWV. 

30 However, the decrease was seen in the data divided into seasonal subsets, and associations 

31 between cloud fraction and aerosol were consistent with expectations that follow from the 

32 aerosol-cloud fraction feedback hypothesized by Koren et al. (2008). At lower cloud 
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1 fractions when the aerosol has greater exposure to sunlight, the decrease of cloud fraction 

2 attributed to aerosol absorption effects increased in magnitude and was shifted to lower values 

3 of AOD. 

4 Relationships between cloud optical depth and aerosol loading are more difficult to interpret. 

5 Plots of COD against AOD showed an initial increase of cloud optical thickness, and then a 

6 turning point where clouds appeared to become optically thinner as aerosols increased. 

7 Increasing COD with aerosol can be attributed physically to the processes described by 

8 Twomey (1977), but the decrease of COD with AOD is best described by a combination of 

9 physical processes and satellite retrieval artifacts. Absorbing aerosol may cause cloud 

10 droplets to evaporate and clouds to thin. This is a legitimate physical process. On the other 

11 hand, the dark aerosol in and above the cloud may also decrease the cloud reflectance 

12 observed by the satellite, which the satellite retrieval interprets as a decrease in COD. As the 

13 cloud thins, subpixel holes in the cloud open, allowing darker reflectance from the surface 

14 beneath to again darken the cloud reflectance measured by the satellite. Again the cloud 

15 reflectance is artificially too low, resulting in COD retrievals that are too low, which 

16 contributes to the strong decrease of COD with AOD observed in the data analysis. 

17 Assuming that the decreasing trend in COD with AOD is due to both radiative effects as well 

18 as retrieval artifacts, our results suggest that GCMs that do not include treatment of aerosol 

19 absorption in clouds are inaccurate over regions with high concentrations of absorbing 

20 aerosol. 

21 This study also addressed the effect of aerosols on CWV and clouds over different land cover 

22 types. Land cover was segregated into forested and deforested groupings using the MODIS 

23 land cover product, and the response of CWV to increasing AOD was analyzed between the 

24 forested and deforested land. We found that aerosols do not have a noticeable effect on 

25 CWV. Thus, it can be assumed that temporal changes in CWV are largely a function of the 

26 synoptic-scale meteorology and that aerosol is completely independent of the column water 

27 vapor at the local level. The independence of CWV and AOD is an inherent assumption in 

28 the former portion of the paper, supported by the analysis in the latter portion of the paper. 

29 This study showed that the relationship between cloud fraction and AOD is quite sensitive to 

30 land cover type. Microphysical effects appeared to be stronger over the deforested land 

3 1 compared to the forested land, increasing cloud fraction with AOD to a greater extent over the 

32 deforested land compared to the forest. The difference in the aerosol microphysical effect 
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1 over forested and deforested land appears to be linked to the differing lower tropospheric 

2 stability over these two surfaces, but a full understanding of the complex interaction between 

3 aerosols, clouds and surface types cannot be achieved from the statistical approach used here. 

4 A 3D cloud resolving model with adequate representation of surface-atmosphere exchange 

5 and ability to similate mesoscale circulations and cloud microphysical processes will be 

6 required. 

7 The relationships shown here cannot, with certainty, be extrapolated to different times of the 

8 year or different regions of the world. These aerosol-cloud correlations are largely dependent 

9 on the physical and optical characteristics of the biomass burning aerosol (i.e. size 

10 distribution, composition); therefore a different aerosol mixture may or may not have a 

11 different impact on clouds. However, a similar analysis to this one over other biomass 

12 burning regions or heavily fossil-fuel polluted regions would be helpful in determining if 

13 these relationships could be extended more generally. More importantly, modeling-based 

14 studies are required to assign causation to the correlations found here. 

15 
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1 Table 1. Average cloud and atmospheric profile statistics for all warm cloud retrievals for 

2 both forested and deforested land in the study region. The averaging period is for August 

3 through October for the years 2004, 2005, 2006, and 2007. 



Forest 

Deforested 

Cloud Effective Radius (um) 

15.9 

14.7 

Cloud Optical Depth (-) 

8.95 

8.07 

Cloud Fraction (-) 

0.52 

0.63 

Cloud Top Pressure (hPa) 

796.4 

817.9 

Cloud Water Path (g/m 2 ) 

81.3 

68.0 

Column Water Vapor (cm) 

3.76 

3.54 

850 hPa Temperature (K) 

292.6 

293.0 

1 000 hPa Temperature (K) 

307.2 

311.4 

Low-Level Stability 

-14.6 

-18.4 

(Temperature 850 hPa minus 



Temperature 1000 hPa (K) ) 



AOD at 550 nm (-) 

0.92 

0.89 
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2 Figure 1. Land cover classifications for the year 2004. Green pixels represent evergreen 

3 forests and yellow pixels represent pasture or deforested land. Blue pixels represent land 

4 classifications that were not included in either of these two categories. The 5° x 5° study 

5 region is outlined with a red box. 
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Date 


1 

2 Figure 2. (a) Column water vapor and cloud fraction binned by day of the year for the months 

3 of August through October, for the years 2004, 2005, 2006, and 2007. Aerosol optical depth 

4 at 550 nm binned by day of the year for the same time period as in (a). Data are accumulated 

5 into eight bins of equal width between August 1 st and October 3 1 st . 
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1 

2 Figure 3. (a) Column water vapor binned by aerosol optical depth for all co-located warm 

3 cloud retrievals for the months of August to October, for the years 2004, 2005, 2006, and 

4 2007. Error bars denote the standard error in each bin. (b) Cloud fraction binned by aerosol 

5 optical depth for all warm cloud retrievals for the same time period as in (a), (c) Column 

6 water vapor binned by aerosol optical depth for all co-located warm cloud retrievals with 

7 cloud fractions less than 0.5. (d) Cloud fraction binned by aerosol optical depth for all warm 

8 cloud retrievals with cloud fractions below 0.5. 
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1 

2 Figure 4. (a) Column water vapor binned by aerosol optical depth for all co-located warm 

3 cloud retrievals at varying percentiles of column water vapor for the months of August to 

4 October, for the years 2004, 2005, 2006, 2007. (b) Cloud fraction binned by aerosol optical 

5 depth for all warm cloud retrievals for the same time period is in (a). 
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1 

2 Figure 5. (a) Cloud optical depth binned by aerosol optical depth for all warm cloud retrievals 

3 at varying percentiles of column water vapor for the months of August to October, for the 

4 years 2004, 2005, 2006 and 2007. (b) Same as (a) but only for retrievals with cloud top 

5 pressures between 800 hPa and 850 hPa. 
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Aerosol Optical Depth (-) Aerosol Optical Depth (-) 


1 

2 Figure 6. (a) Difference between the average sampled day of the year in each aerosol bin for 

3 the pasture and forest sites. The sampling algorithm ensures that the same number of samples 

4 per day is used from each land cover type in each AOD bin. In panels a, b and c, red lines 

5 represent the results of five random simulations, where forest and pasture pixels are mixed. 

6 Data represent the months August to October, for the years 2004 to 2007. (b) Differenced 

7 column water vapor between pasture and forest land cover types binned by AOD. (c) 

8 Differenced cloud fraction between pasture and forest land cover types binned by AOD. (d) 

9 Black line: same as (c), red and blue lines: Stratification by low level stability, defined as the 

10 temperature at 850 hPa minus the temperature at 1000 hPa. The stratification values are 

1 1 different between pasture and forest for each stratification case. 
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a) 



Aerosol Optical Depth AERONET 500 nm (-) 



1 

2 Figure SI. (a) Comparison of co-located MODIS AOD retrievals with automatic sun/sky 

3 radiometers of the Aerosol Robotic Network (AERONET) at Abracos Hill (62.358° W, 

4 10.760° S) and Ji Parana SE (61.852° W, 10.934° S) for days used in the study between the 

5 years 2004 through 2007 during the months of August through October. These two sites are 

6 within our 5° x 5° study region, (b) Comparison of co-located MODIS CWV retrievals with 

7 automatic AERONET retrievals for the same time period as in (a). Correlation coefficients 

8 between MODIS and AERONET data are also included in each plot. 


38 


Missing/Bad 
No Signal 
Sub-Surface 
Surface 
Stratosphere 
Aerosol 
Cloud 
Clean Air 

(-60.7,-13.0) (-60.9,-12.0) (-61.1,-11.0) (-61.3,-10.0) (-61.6, -9.0) (-61.8, -8.0) 

Location (Lon.Lat) 



□u 



• 

0 

Uncertainty < 1 0% 

45 



• 

o 

Uncertainty < 20% 



* s 

• 

Uncertainty < 30% 

40 



* * * * 

• 

Other 


S' 35 

Q. 

Q 30 
S 25 

Cl 

° 20 
"O 

10 

5 

0 


°o°o <P 

op 

9 <V 
o ft 


o 

V o 

8 o 

% o e ° 


o° 

o ° 
o 


°o<?° 


* •*> 
o OO u o 

oo g o ° 8 ° 
„ 0 < 


Vo 


; ooo 

(D°°c 


o 

0 


ft) CO 

ft o 


o 

o o 


cP 


O o 
Po 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Aerosol Optical Depth (-) 


0.8 0.9 


1 

2 Figure S2. (a) Vertical Feature Mask from the CALIPSO lidar on a path through the study 

3 region for a single day during the study period, August 12 , 2006. On this day, clouds mostly 

4 form near the top of the aerosol layer. Clouds that extend throughout the column may be 

5 misclassified heavy aerosol plumes, (b) Scatter plot of cloud optical depth versus aerosol 

6 optical depth on the same day as (a), throughout the 5° x 5° study region. Colors correspond 

7 to the cloud optical depth uncertainty parameter, defined as a percentage of the retrieved 

8 value. 
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